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Abstract. An efficient algorithm is presented to compute the character¬ 
istic polynomial of a threshold graph. Threshold graphs were introduced 
by Chvatal and Hammer, as well as by Henderson and Zalcstein in 1977. 
A threshold graph is obtained from a one vertex graph by repeatedly 
adding either an isolated vertex or a dominating vertex, which is a ver¬ 
tex adjacent to all the other vertices. Threshold graphs are special kinds 
of cographs, which themselves are special kinds of graphs of clique-width 
2. We obtain a running time of 0(n log 2 n ) for computing the character¬ 
istic polynomial, while the previously fastest algorithm ran in quadratic 
time. 
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1 Introduction 

The characteristic polynomial of a graph G = (V, E ) is defined as the char¬ 
acteristic polynomial of its adjacency matrix A, i.e. x(G, A) = det(A I—A). 
The characteristic polynomial is a graph invariant, i.e., it does not depend 
on the enumeration of the vertices of G. The complexity of computing the 
characteristic polynomial of a matrix is the same as that of matrix mul¬ 
tiplication [9,12] (see [1, Chap.16]), currently 0(n 2 ' 3 ' 6 ) [3]. For special 
classes of graphs, we expect to find faster algorithms for the character¬ 
istic polynomial. Indeed, for trees, a chain of improvements [15,11] re¬ 
sulted in an 0(n log 2 n) time algorithm [6]. The determinant and rank of 
the adjacency matrix of a tree can even be computed in linear time [4], 
For Threshold graphs (defined below), Jacobs et al. [8] have designed an 
0(n 2 ) time algorithm to compute the characteristic polynomial. Here, we 
improve the running time to 0(n log 2 n). As usual, we use the algebraic 
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complexity measure, where every arithmetic operation counts as one step. 
Throughout this paper, n = \V\ is the number of vertices of G. 

Threshold graphs [2,7] are defined as follows. Given n and a sequence 
b = 6 n _i) G {0, the threshold graph Gf, = (V, E) is defined 

by V = [ra] = {1,..., n}, and for all i < j, {i,j} G E iff bi = 1. Thus Gb 
is constructed by an iterative process starting with the initially isolated 
vertex n. In step j > 1, vertex n — j + 1 is added. At this time, vertex 
j is isolated if bj is 0 , and vertex j is adjacent to all other (already 
constructed) vertices {j + 1,..., n} if bj = 1. It follows immediately that 
Gb is isomorphic to Gy iff b = b'. Gb is connected if b\ = 1 , otherwise 
vertex 1 is isolated. Usually, the order of the vertices being added is 
1,2,... ,n instead of n, n — 1,..., 1. We choose this unconventional order 
to simplify our main algorithm. 

Threshold graphs have been widely studied and have several applica¬ 
tions from combinatorics to computer science and psychology [ 10 ]. 

In the next section, we study determinants of weighted threshold 
graph matrices, a class of matrices containing adjacency matrices of thresh¬ 
old graphs. In Section 3, we design the efficient algorithm to compute the 
characteristic polynomial of threshold graphs. We also look at its bit com¬ 
plexity in Section 4, and finish with open problems. 


2 The determinant of a weighted threshold graph matrix 


We are concerned with adjacency matrices of threshold graphs, but we 
consider a slightly more general class of matrices. We call them weighted 
threshold graph matrices. Let '"^ n be the matrix with the following 

entries. 


m: 


did2---d n 


bib2---bn-i 


bi if i < j 

< bj if j < i 

y di if i = j 


Thus, the weighted threshold matrix for (6162 • • • b n - 1 ; d\d 2 ■ ■ ■ d r 
like this. 


looks 


l\ / rdid2...d n 


U,J[ VL Ui . 

b\ d 2 b 2 . 

.. b 2 

b 2 

b\ b 2 d 3 . 

• • b 3 

b 3 

bi b 2 63 . 

• • d n —\ b n — 1 

\h b 2 63 . 

• • b n — 1 d n J 
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In order to compute the determinant of M d ^ d2 "' dn , we subtract the 
penultimate row from the last row and the penultimate column from 
the last column. In other words, we do a similarity transform with the 
following regular matrix 


i.e., 


/1 0 0 ... 

0 

0 ^ 

0 10... 

0 

0 

00 1... 

0 

0 

00 0... 

1 

-1 

^00 0... 

0 

1 / 

if i = j 




P — 


— 1 if i 


n and j = n — 1 


0 otherwise. 


The row and column operations applied to iV /^ 1 ^ 2 '" dn i produce the similar 
matrix 


P T M?\ d2 "; dn P = 

0lb2-..0n-l 


/d\ b\ b\ .. 
b\ d 2 b 2 .. 
bi b 2 d 3 .. 


bi 

b 2 

b 3 


0 \ 
0 

0 


b\ b 2 b 3 . . . d n ~i b n —i d n —\ 

\ 0 0 0 ... b n — i d n — i d n T d n —\ 2 b n _\J 


Naturally, the determinant of P is 1, implying 


det (P T P) = del («, 


di d 2 ..-dn 

6i6 2 ...6 n _i 


Furthermore, we observe that P T M^ 1 ^ 2 "^ n P has a very nice pattern. 


P T M^ 2 "i dn P = 

Oib2...b n -i 

( 

l\/[dl d2 • • • d n — 1 
iV %6 2 ...6 n _ 2 

bi \ 

b 2 

b 3 




bji—i d n —1 


\ 62 ^3 • • • bn—i d n —i d n -b d n —i 2 b n -i) 
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To further compute the determinant of P T P, we use Laplacian 

expansion by minors applied to the last row. 


** (<■£::£:,) = de * { pT Ktt-, p ) 

= (d n + <i„_, - 2f» n _i) det - d*.,) 2 det 


By defining the determinant of the 0x0 matrix with n = 0 

to be 1, and checking the determinants for n = 1 and n = 2 directly, we 
obtain the following result. 

Theorem 1. = det {^^2 b n 1) determined by the recurrence 

equation 


D 


n 


1 if n = 0 

< di if n = 1 

k ( dn T d n —1 < 2b n —\)D n —\ (6 n -i d n —\)~D n —2 if n. 2 


□ 

This has an immediate implication, as we assume every arithmetic oper¬ 
ation takes only 1 step. 

Corollary 1. The determinant of an n x n weighted threshold graph ma¬ 
trix can be computed in time 0(n). 

Proof. Every step of the recurrence takes a constant number of arithmetic 
operations. 

□ 

For arbitrary matrices, the tasks of computing matrix products, ma¬ 
trix inverses, and determinants are all equivalent [1, Chap. 16], currently 
0(n 2,376 ) [3]. For weighted threshold graph matrices, they all seem to 
be different. We have just seen that the determinant can be computed 
in linear time, which is optimal, as this time is already needed to read 
the input. The same lower bound holds for computing the characteristic 
polynomial, and we will show an 0(n log 2 n) algorithm. It is not hard to 
see that the multiplication of weighted threshold graph matrices can be 
done in quadratic time. This is again optimal, because the product is no 
longer a threshold graph matrix, and its output requires quadratic time. 
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3 Computation of the Characteristic Polynomial of a 
Threshold Graph 

The adjacency matrix A of the n- vertex threshold graph G defined by the 
sequence (b \,..., b n -\) is the matrix _ i5 and the characteristic 

polynomial of this threshold graph is 

x(G, A) = det(AI - A) = det . 

This immediately implies that any value of the characteristic polynomial 
can be computed in linear time. 

The characteristic polynomial itself can be computed by the recur¬ 
rence equation of Theorem 1. Here all di = A, and D n , as the characteris¬ 
tic polynomial of an n-vertex graph, obviously is a polynomial of degree 
n in A. Now, the computation of D n from ZA„_i and D n _ 2 according to 
the recurrence equation is a multiplication of polynomials. It takes time 
O(n), as one factor is always of constant degree. The resulting total time 
is quadratic. The same quadratic time is achieved, when we compute the 
characteristic polynomial x(G, A) for n different values of A and interpo¬ 
late to obtain the polynomial y(G, A). 

We want to do better. Therefore, we write the recurrence equation of 
Theorem 1 in matrix form. 

f D n A / d n T- d n —1 26 n _i (b n —1 — ( Z) n _iA 

\D n -i) V 1 0 J \Dn_2j 

Noticing that Dq = 1 and D\ = A, and all di = A, we obtain the following 
matrix recurrence immediately. 

Theorem 2. For 

Bn = ^ A-M -(^-A) 2 A 

we have 

□ 

This results in a much faster way to compute the characteristic polynomial 
X(G, A). 

Corollary 2. The characteristic polynomial y(G, A) of a threshold graph 
G with n vertices can be computed in time 0(nlog 2 n). 
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Proof. For every i, all the entries in the 2x2 matrix Bi are polynomials 
in A of degree at most 2. Therefore, products of any k such factors have 
entries which are polynomials of degree at most 2k. To be more precise, 
actually the degree bound is k, because by induction on k, one can easily 
see that the degrees of the i, j-entry of such a matrix is at most 

k for i = 1 and j = 1, 
k + 1 for i = 1 and j = 2, 
k — 1 for i = 2 and j = 1, 
k for i = 2 and j = 2, 


But the bound of 2k is sufficient for our purposes. W.l.o.g., we may assume 
that n — 1 (the number of factors) is a power of 2. Otherwise, we could 
fill up with unit matrices. Now the product B n _\B n _2 - - - -Si is computed 
in log(n — 1) rounds of pairwise multiplication to reduce the number of 
factors by two each time. In the rth round, we have n2~ r pairs of matrices 
with entries of degree at most 2 r , requiring 0(n2~ r ) multiplications of 
polynomials of degree at most 2 r . With FFT (Fast Fourier Transform) 
this can be done in time 0(nr). Summing over all rounds r results in a 
running time of 0(n log 2 n). 

□ 

Omitting the simplification of di = A in Theorem 2, we see imme¬ 
diately, that also the characteristic polynomial of a weighted threshold 
graph matrix can be computed in the same asymptotic time of 0(n log 2 n). 

4 Complexity in the Bit Model 

By definition, the characteristic polynomial of an n-vertex graph can be 
viewed as a sum of n! monomials with coefficients form {—1,0,1}. Thus 
all coefficients of the characteristic polynomial have absolute value at 
most nl, and can therefore be represented by binary numbers of length 
O(relogn). The coefficients can indeed be so big. An example is the con¬ 
stant term in the characteristic polynomial of the clique K n . Its absolute 
value is the number of derangements (permutations without fixed points), 
which asymptotically converges to n!/e. 

With such long coefficients, the usual assumption of arithmetic oper¬ 
ations in linear time is actually unrealistic for large n. Therefore, the bit 
model might be more useful. We can use the Turing machine time, be¬ 
cause our algorithm is sufficiently uniform. No Boolean circuit is known 
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to compute such things with asymptotically fewer operations than the 
number of steps of a Turing machine. 

We can use the fast m log m2°( log integer multiplication algorithm [5] 
(where m is the length of the factors) to compute the FFT for the poly¬ 
nomials. A direct implementation, just using fast integer multiplication 
everywhere, results in time 


(nr)(2V 2 2°( log ‘ r )) = n2 r r 3 2°( log ‘ r ) 


for the rth round where 0(n2~ r ) pairs of polynomials of degree 0(2 r ) are 
multiplied. The coefficients of these polynomials have length 0(2 r r). As 
the coefficients and the degrees of the polynomials increase at least geo¬ 
metrically, only the last round with r = log n counts asymptotically. The 
resulting time bound is 0(n 2 log 3 n2°( log n )). Using Schonhage’s [13] idea 
of encoding numerical polynomials into integers in order to do polynomial 
multiplication, a speed-up is possible. Again only the last round matters. 
Here a constant number of polynomials of degree 0{n) with coefficients 
of length O(nlogn) are multiplied. For this purpose, each polynomial is 
encoded into a number of length 0(jn? logn), resulting in a computation 
time of 

n 2 log 2 n 2°( log * n ). 

Actually, because the lengths of coefficients are not smaller than the 
degree of the polynomials, no encoding of polynomials into numbers is 
required for this speed-up. In this case, one can do the polynomial mul¬ 
tiplication in a polynomial ring over Fermat numbers as in Schonhage 
and Strassen [14]. Then, during the Fourier transforms all multiplications 
are just shifts. Fast integer multiplication is only used for the multipli¬ 
cation of values. This results in the same asymptotic n 2 log 2 n 2°( log n ) 
computation time with a better constant factor. 

5 Open Problems 

We have improved the time to compute the characteristic polynomial of a 
threshold graph from quadratic to almost linear (in the algebraic model). 
The question remains whether another factor of logn can be removed. 
More interesting is the question whether similarly efficient algorithms 
are possible for richer classes of graphs. Of particular interest are larger 
classes of graphs containing the threshold graphs, like cographs, graphs of 
clique-width 2, graphs of bounded clique-width, or even perfect graphs. 
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